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Summary 

The Plasmodium falciparum genome is equipped with 
several subtelomeric gene families that are implicated 
in parasite virulence and immune evasion. Members of 
these families are uniformly positioned within hetero- 
chromatic domains and are thus subject to variegated 
expression. The best-studied example is that of the 
var family encoding the major parasite virulence 
factor P. falciparum erythrocyte membrane protein 1 
(PfEIUIPI). PfEMPI undergoes antigenic variation 
through switches in mutually exclusive i^argene tran- 
scription. 1/arpromoters function as crucial regulatory 
elements in the underlying epigenetic control strategy. 
Here, we analysed promoters of upsA, upsB and upsC 
var, rifA 7-type rif, stevor, phist and pfmc-2tm genes 
and investigated their role in endogenous gene tran- 
scription by comparative genome-wide expression 
profiling of transgenic parasite lines. We find that the 
three major i/ar promoter types are functionally equal 
and play an essential role in singular gene choice. 
Unlike i/ar promoters, promoters of non-var families 
are not silenced by default, and transcription of non- 
i/ar families is not subject to the same mode of mutu- 
ally exclusive transcription as has been observed for 
i^ar genes. Our findings identified a differential logic in 
the regulation of varand other subtelomeric virulence 
gene families, which will have important implications 
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for our understanding and future analyses of pheno- 
typic variation in malaria parasites. 

Introduction 

The blood stages of Plasmodium falciparum evade 
antibody-mediated host immunity by altering parasite- 
encoded antigens exposed on the surface of infected red 
blood cells (iRBCs) through antigenic variation (Biggs 
etal., 1991; Roberts etal., 1992; Smith etal., 1995; Kyes 
et al., 2001 ). Antigenic variation results from the amplifica- 
tion of hypervarlable gene families, coupled to the ability to 
switch expression of Individual genes. Notably, these gene 
families are predominantly located near the telomeres 
(Gardner et al., 2002). The positional clustering of antigen 
gene families in subtelomeric regions Is not unique to 
P. falciparum but Is common to a wide range of pathogenic 
protozoa and fungi (Myler etal., 1984; Barry etal., 2003; 
De Las etal., 2003; Keely etal., 2005; Deltsch etal., 
2009). Subtelomeric gene families show high rates of 
recombination to facilitate diversity and this has been 
demonstrated In the model yeast Schizosaccharomyces 
pombe (Cooper etal., 1998) and In P. falciparum, where 
meiotic recombination among var genes is at least eight 
times more frequent than the estimated genomic average 
(Freltas-Junlor etal., 2000; Taylor etal., 2000). Proximity 
to telomeres also modulates gene expression as exempli- 
fied by variegated silencing of genes nearby telomeres 
(Gottschling etal., 1990). Indeed, P. falciparum subte\orr\- 
eric regions are heterochromatic and this repressive envi- 
ronment facilitates silencing and variegated gene 
expression (Flueck etal., 2009; Lopez-Rubio etal., 2009; 
Salcedo-Amaya etal., 2009). 

Subtelomeric P. falciparum gene families include var 
(Baruch etal., 1995; Smith etal., 1995; Su etal., 1995), rif 
(repetitive Interspersed family) (Weber, 1988), stevor{sub- 
telomerlc variable open reading frame) (Cheng etal., 
1998), phist {Plasmodium helical Interspersed subtelom- 
eric) (Sargeant et al., 2006), pfmc-2tm (P. falciparum Mau- 
rer's clefts two transmembrane) (Sam-Yellowe et al. , 2004) 
and the hyp families 1-17 (Sargeant etal., 2006). The 
best-studied example Is that of the ~ 60 member var gene 
family, var gene transcription Is controlled by conserved 
promoter sequences that can be subgrouped into three 
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main classes according to sequence and cliromosomal 
location (Voss ef a/., 2000; Gardner etal, 2002; Kraemer 
and Smith, 2003; Lavstsen ef a/., 2003). PfEMPI , encoded 
by i^ar genes, is exposed on the iRBC surface and medi- 
ates adherence to various endothelial receptors in the 
human host (Leech ef a/., 1 984; Su ef a/., 1 995; Kyes ef a/., 
2001). The resulting sequestration of RBC aggregates in 
the microvasculature is a major cause of severe clinical 
outcomes such as cerebral and pregnancy-associated 
malaria (MacPherson efa/., 1985; Carlson efa/., 1990; 
Pongponratn efa/., 1991; Rowe efa/., 1995). Moreover, 
antigenic variation of PfEMPI contributes substantially to 
chronic disease and transmission. 

Unlike PfEMPI, the function of other exported protein 
families in parasite biology and disease is unclear, but 
there is indication for roles in host-parasite interactions 
(Fernandez efa/., 1999; Kyes efa/., 1999; Abdel-Latif 
efa/., 2002; 2003; Lavazec efa/., 2007; Schreiber efa/., 
2008). With over 130 members, the r/f gene family is the 
largest in P. falciparum (Kyes efa/., 1999; Joannin efa/, 

2008) . A-type RIFINs are associated with parasite-induced 
membranous structures in the iRBC cytosol termed 
Maurer's clefts (MC) and the infected RBC surface. B-type 
RIFINs have been shown to reside within the parasite 
(Petter efa/., 2007; 2008). Members of the ~ 35-member 
STEVOR family localize to the surface of infected erythro- 
cytes and are also associated with MCs (Lavazec efa/., 
2006; BIythe ef a/., 2008; Khattab ef a/., 2008; Niang ef a/., 

2009) . The PFMC-2TM family (12 members) are associ- 
ated with MCs and the RBC surface (Sam-Yellowe efa/., 
2004; Lavazec efa/., 2006). PRIST and HYP proteins are 
predicted to be exported to the IRBC (Sargeant efa/, 
2006), but this has not yet been experimentally confirmed. 

var genes are transcribed in a mutually exclusive 
manner such that in single parasites only one uar gene is 
active and all other members are silenced (Scherf efa/., 
1998). However, simultaneous transcription of two var\oc\ 
has been observed in wild-type and transgenic parasites 
(Dzikowski efa/., 2007; Brolin efa/., 2009; Joergensen 
efa/., 2010). var promoters are essential in establishing 
alternative states of activity (Dzikowski efa/., 2006; 2007; 
Voss efa/., 2006; 2007; Chookajorn efa/., 2007; Lopez- 
Rubio ef a/., 2007). Transfection experiments using episo- 
mal promoters to drive expression of drug selectable 
markers proved a great strategy to investigate transcrip- 
tional regulation of i/argenes (Dzikowski ef a/., 2006; Voss 
efa/., 2006; 2007). Upon transfection, var promoters are 
silenced by default (Calderwood efa/., 2003; Frank efa/, 
2006; Voss efa/., 2006; 2007) and this silenced state 
appeared to require the presence of the var gene intron 
(Deitsch efa/., 2001; Frank efa/, 2006). Subsequent 
studies showed that the intron is not strictly required for 
silencing; rather the c/s pairing of a var promoter with the 
intron'sown promoter activity, or that of other heterologous 



promoters, is important (Dzikowski efa/., 2007). Interest- 
ingly, drug-induced selection for active promoters is suffi- 
cient to infiltrate drug-selectable markers into the pathway 
of mutually exclusive transcription (Dzikowski efa/., 2006; 
Voss efa/., 2006; 2007). In such parasites, the active 
episomal promoters cause a shutdown of endogenous var 
transcription and, consequently, absence of PfEMPI on 
the iRBC surface. These studies have been instrumental 
not only in understanding the regulatory mechanisms 
underlying mutually exclusive transcription, but also in 
reconstructing earlier observations with regard to the vari- 
able phenotype of IRBCs (Voss efa/., 2006; D'Ombrain 
efa/, 2007; Elliott efa/., 2007; Maier efa/., 2008). 

Despite most other multigene families in P. falciparum 
being known for many years, our understanding of the 
processes that regulate their transcription is limited. 
Expression of r/f, stevor and pfmc-2tm genes is clonally 
variant In vitro and In vivo (Fernandez ef a/, 1 999; Lavazec 
efa/., 2007; Petter efa/., 2007; 2008; BIythe efa/., 2008; 
Mackinnon ef a/., 2009). Analysis of parasites selected for 
monogenic var expression provided no evidence for 
co-regulation between the active var and neighbouring rif 
or sfevor genes (Sharp efa/., 2006; Cabral and Wunder- 
lich, 2009). These studies indicate that juxtaposition to 
active var genes does not predispose the activity of other 
gene family members. However, upsAand rlfA1 genes, a 
subset of r/f genes positioned in head-to-head orientation 
with upsA var genes, were both upregulated in parasites 
selected to express PfEMPI variants associated with 
severe disease in children (Wang ef a/., 2009). Moreover, a 
recent study observed interference of several episomal 
gene family promoters with transcription of endogenous 
gene families suggesting that multiple gene families may 
share a common transcriptional control mechanism (Howitt 
efa/., 2009). 

Here, we used genome-wide comparative transcrip- 
tional profiling to investigate the extent to which activation 
of extrachromosomal virulence gene promoters affects 
transcription of endogenous gene families. Our results 
establish a clear functional distinction between var and 
non-varpromoters. They further demonstrate that varpro- 
moters of all types play an essential part in singular gene 
choice, and that transcription of non-vargene families is 
not subject to the same mode of transcriptional regulation 
that controls mutually exclusive transcription of vargenes. 

Results 

Generation of stable reporter parasite lines 

Parasites of the strain 3D7 were transfected with con- 
structs carrying eight different promoters: upsA-type var 
(PF13_0003), upsB-type var (PFL0005w), upsC-type var 
(PFL1960W), rlfA1 (PF13_0004), stevor (PFL2610w), 
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phistb (PFL2540w), pfmc-2tm (final gene ID unknown due 
to high sequence conservation between pfmc-2tm promot- 
ers) and cam (PF14_0323) (control). Note that the rifA1 
and upsA promoters are present on the sense and anti- 
sense strand, respectively, of the 2.8 kb intergenic 
sequence separating these head-to-head oriented genes. 
All sequences were cloned into the context of the parental 
reporter vector pBcam by replacement of the cam pro- 
moter (Fig. 1A). The blasticidin-S deaminase {bsd) resis- 
tance cassette was used to obtain transfectants carrying 
stable episomes. Importantly, this initial selection is inde- 
pendent of the activity of the test promoters driving tran- 
scription of the downstream drug-selectable reporter hdhfr- 
gfp (human dihydrofolate reductase fused to green 
fluorescent protein). All constructs carry a var gene intron 
element downstream of the hdhfr-gfp cassette for consis- 
tency. We chose this approach, rather than cloning family- 
specific introns, since the role of the var intron in silencing 
and mutual exclusion is exerted by the intron's own pro- 
moter activity (i.e. a functional feature of the intron) and is 
thus sequence-independent (Dzikowski etal., 2007). 

The BSD-selected reporter lines 3D7/upsAt, 3D7/upsB, 
3D7/upsC, 3D7/rifA1, 3D7/stevor, 3D7/phist, 3D7/pfmc- 
2tm and the control line 3D7/cam exhibited similar plasmid 
copy numbers ranging from two to eight copies per parasite 
(Fig. 1 B). After WR9921 0 (WR) selection, average plasmid 
numbers increased with the highest number of 33 
observed in 3D7/rifA1. Unexpectedly, 3D7/upsAt was 
completely refractory to WR selection in seven challenge 
experiments on three independently generated BSD- 
resistant populations. The reason for failure to select for 
active upsA promoters is unknown but may have been 
related to efficient upsA silencing and/or insufficient provi- 
sion of stable hc//7frtranscripts to confer WR resistance. In 
an attempt to circumvent these two possible obstacles, we 
exploited the possibility that rifA1 promoter activity may 
predispose the head-to-head upsA promoter for activation. 
We therefore replaced the entire bsd cassette in pBrifAl 
with an inverted bsd gene and terminator element to create 
pBupsA (Fig. 1C). In this context, the reverse complement 
sequence of the rifA1 upstream region acts as an upsA 
promoter driving fesd transcription. Indeed, parasites trans- 
fected with this construct were successfully selected on 
BSD. Hence, this BSD-resistant population (3D7/upsA) 
carried an activated upsA promoter and was used in all 
subsequent experiments. Noteworthy, populations trans- 
fected with pBupsA were also successfully selected on 
WR, and reciprocal drug swaps had no effect on parasite 
growth suggesting that both activities of this bi-directional 
promoter were active simultaneously (data not shown). 
Southern blot analysis revealed that plasmids were main- 
tained episomally in all lines except 3D7/stevor where 
integration into the endogenous stevor locus occurred 
(Fig.SI). 
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Fig. 1. Transgenic reporter cell lines. 

A. Schematic map of the parental plasmid pBcam. The bsd 
resistance cassette (white box) selects for stable plasmid 
maintenance. The hdhfr-gfp cassette (grey box) allows selecting for 
active promoters of interest (red arrow). The var gene Intron is 
indicated by a v-shape. A 0.5 kb TARE6 element is shown by 
vertical lines. hsp86 5', hsp86 promoter; Pb DT 3', P. berghei 
d/ifr-thymldilate synthase terminator; firp2 3', histidlne-rich protein 2 
terminator 

B. Descendants of pBcam were obtained by replacing the cam 
promoter with promoters of interest (red arrows) using Bglll and 
Notl. Numbers represent nucleotide positions in relation to the ATG 
start codon. Grey dashed arrows represent the bi-directional 
upsA/r/f/A? promoter. Plasmid names and cell lines are indicated. 
White and grey bars represent average plasmid copy numbers in 
unselected and WR-selected cell lines respectively, n.a., not 
available. 

C. Top: schematic map of pBupsA where the bsd and hdhfr-gfp 
selectable markers are controlled by the bi-directional upsA/riM) 
promoter. Bottom: average plasmid copy number in BSD-selected 
3D7/upsA parasites. 
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Initial assessment of overall transcript profiles at four 
consecutive intervals during intra-erythrocytic 
development 

We harvested total RNA at four consecutive time points 
[TP1: 6-14 11 post invasion (hpi); TP2: 14-22 lipi; TPS: 
22-30 Up\; TP4: 30-38 hpi] during the intra-erythrocytic 
developmental cycle (IDC) from 3D7 wild-type parasites 
and each of the WR-selected transfected lines (except for 
3D7/upsA where a BSD-selected population was used). 
This sampling strategy was designed to Include the tem- 
poral peaks of transcription for each endogenous gene 
family (Scherf etal., 1998; Kyes etai., 1999; 2000; Kavi- 
ratne et aL, 2002; Sam-Yellowe et aL, 2004). Relative tran- 
script levels for each gene were determined relative to a 
3D7 reference cDNA pool by hybridization to a genome- 
wide long oligonucleotide microarray (Hu etal., 2007) 
(GEO accession GSE31829). Hierarchical average 
linkage clustering revealed four distinct clusters corre- 
sponding to the four sampling time points (Fig. 2A). Within 
each time point, the samples showed a high degree of 
correlation but the individual datasets clustered differently. 
The overall similarity between individual transcriptomes 
was lowest in early ring stages and increased gradually 
with progression through the IDC (Fig. 2B). This is 
explained by the variegated expression of a large number 
of genes expressed specifically in ring-stage parasites 
(Bozdech etal., 2003; IVIarti etal., 2005; Llinas etal, 
2006). Compared with both control parasites (3D7 wild 
type and 3D7/cam), the ring-stage transcriptomes of 3D7/ 
upsA, 3D7/upsB, 3D7/upsC and 3D7/rifA1 displayed a 
higher number of deregulated genes than 3D7/stevor, 
3D7/phist and 3D7/pfmc-2tm. Notably, all reporter line 
transcriptomes consistently showed a higher degree of 
similarity to the control transfectant 3D7/cam than to 
untransfected parasites (Fig. 2B). We explain this impor- 
tant observation by the likelihood that perturbations asso- 
ciated with the generation of transgenic parasites per se 
cause alterations in mRNA profiles. Indeed, we found 127 
genes either up- or downregulated by more than twofold in 
at least one time point in all transfected lines compared 
with the untransfected control (Fig. S2). Interestingly, this 
set contained several genes coding for proteins involved in 
DNA/RNA metabolism, protein translation and export, and 
stress response. These changes are compatible with pos- 
sible responses to plasmid maintenance and replication as 
well as to the mode of BSD action, which inhibits protein 
synthesis (Yamaguchi etal., 1965). 

Promoters of var and non-var subtelomeric gene 
families display distinct functional characteristics 

As an important consideration for comparative analyses, 
we tested whether the stage-specific activities of the epi- 
somal promoters matched the temporal transcription 



B 



# genes >4-fotd change # genes >4-fold change 
relative to 3D7 relative to 3D7/cam 



in ■.- 




^g^D7/upsB 

'— „ 3D7/stevor 
097 

|^g^3D7/phist 

— 3D7/pfmc-2tm 
0,97 

^„ 3D7/upsC 
0.96 

'— 3D7 
0,94 

— 3D7/cam 

1^96^°™'*^ 
' — 3D7/upsA 



^9g3D7/nfA1 



-^3D7/phist 

— 3D7/cam 
0,97 

-g3D7/upsC 
^gg3D7/upsA 

' — 3D7 

095 

' — 3D7/pfmc-2tm 
0,93 

— TT 3D7/stevor 
0 93 

3D7/upsB 



P 



[=] 



[=] 



[] 



[] 



Z] 



□ upregulated 

□ downregulated 



[] 



IZI 

c 

LI 



^1 



CI 

u 



CI 



[] 



□ upregulated 

□ downregulated 



Fig. 2. Hierarchical clustering of 36 transcriptomes. 

A. Hierarchical average linkage clustering of 3D7 and eight 
transgenic cell lines at four time points during the IDC. Numbers 
indicate similarity scores of two joined elements by uncentred 
Pearson correlation. TP, time point (TP1, 6-14 hpi; TP2, 14-22 hpi; 
TPS, 22-30 hpi; TP4, 30-38 hpi). 

B. Bars indicate the number of genes displaying more than fourfold 
change in relative expression (open bars, downregulated; grey 
bars, upregulated) in each cell line and time point relative to 3D7 
wild type (left panel) or the 3D7/cam control transfectant (right 
panel). Parasites transfected with var promoter constructs are 
shaded grey. 



profile of the cognate endogenous gene families. To do 
this, we compared reporter gene mRNA abundance pro- 
files in selected lines to the average relative expression of 
the corresponding gene family members (Fig. S3). This 
analysis revealed high correlations for the var, stevor, 
pfmc-2tm and cam promoters demonstrating that the epi- 
somal sequences contained sufficient information to reca- 
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pitulate the temporal expression of endogenous gene 
families. Intrlgulngly, tlie rifA1 promoter correlated with 
endogenous \/ar rather than r/nranscription implying that 
this bi-directional r/f promoter displays a var-specific acti- 
vation pattern. This interesting observation may indicate 
differential regulation of this particular subset of r/f genes 
compared with other members of the family. For the phistb 
promoter we observed divergence in early schizonts where 
episomal promoter activity increased and endogenous 
ph;sf transcription decreased. The reason for this discrep- 
ancy is unclear but may be due to lack of the full comple- 
ment of regulatory c/s-acting sequences on the 2.5 kb 
upstream region cloned into pBphistb, or different mRNA 
stability properties of the reporter compared with endog- 
enous p/7/snranscripts. 

Several studies reported that upsB and upsC promoters 
are regulated by epigenetic mechanisms and adopt either 
a silenced (default) or an active (activated) state 
(Dzikowski etal., 2006; Frank etal., 2006; Voss etal., 
2006; 2007). By quantifying steady-state mRNA levels in 
unselected (default) and drug-selected (activated) popula- 
tions we asked if the promoters of upsA varand subtelom- 
eric non-i/ar families show the same behaviour. The cam 
promoter that is not epigenetically regulated and the upsB 
and upsC varpromoters served as important negative and 
positive controls for silencing respectively. As expected, 
quantitative reverse transcription PGR (qRT-PCR) analy- 
sis revealed high levels of induction for the upsB promoter 
(18.2-fold ± 1.0 s.d.) and upsC promoter (9.4-fold ± 
3.7 s.d.) relative to the negative control promoter cam 
(2.1 -fold ± 0.4 s.d.) (Figs 3 and S4). Both the upsA 
(3.8-fold ± 0.6 s.d.) and rifA1 promoters (7.1-fold ± 
3.4 s.d.; not significant) displayed weaker but still elevated 
levels of induction compared with the control, which may 
be indicative for silencing of both activities of this 
bi-directional promoter in a considerable fraction of the 
unselected populations. In contrast, we found no evidence 
for alternative states of activity for the stevor, phistb and 
pfmc-2tm promoters. These results therefore suggest for 
the first time that promoters of subtelomeric non-i/arfami- 
lies may be unable to recruit transcriptional silencing 
autonomously. This is in distinction to var promoters that 
are characterized by alternative states of activity and 
reversible responsiveness to epigenetic control mecha- 
nisms (Deitsch etal., 2001; Duraisingh etal., 2005; 
Freitas-Junior etal., 2005; Voss etal., 2006; 2007; 
Chookajorn etal., 2007; Dzikowski etal., 2007; Lopez- 
Rubio etal., 2007). 

The capability to Interfere with transcription of 
endogenous gene families is restricted to var 
gene promoters 

Activation of upsB and upsC promoters leads to down- 
regulation of endogenous var transcription (Dzikowski 
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Fig. 3. Gene family promoters other than I'arare not silenced by 
default. Relative reporter transcript levels were determined by 
qRT-PCR [normalized against transcription of arginyl-tRNA 
synthetase (PFL0900c) and adjusted for plasmid copy numbers]. 
The bars represent the average increase in relative promoter 
activity in drug-selected compared with unselected parasites at the 
time of peak transcription (fold increase). Values derive from three 
independent biological replicates (mean ± s.d.) (see also Fig. S4). 
Note that the value for the upsA promoter was calculated by 
dividing the relative bsd transcript levels in 3D7/upsA by the relative 
hdhfr-gfp levels in unselected 3D7/upsAt. The horizontal dashed 
line indicates the level of induction observed for the control 
promoter cam. Significant levels of induction compared with the 
negative control line 3D7/cam are indicated by asterisks (P< 0.05). 
TP, time point. 

etal., 2006; 2007; Voss etal., 2006; 2007). Here, we 
studied the specificity of this cross-talk by comparative 
analysis of relative i/ar transcription levels in all cell lines. 
As shown in Fig. 4A, mean relative var expression was 
strikingly reduced in early and late ring stages in 3D7/ 
upsA, 3D7/upsB and 3D7/upsC compared with all other 
lines. Interestingly, this was also evident in 3D7/rifA1 to a 
level comparable to that induced by active varpromoters. 
As expected, var transcription was repressed/absent in 
trophozoites and early schizonts (TP3 and TP4) in all lines 
investigated, underscoring the specificity of the changes 
observed for TP1 and TP2. 

The box plot in Fig. 4A compares the distribution of 
relative var expression levels at TP1 . With the exception 
of 3D7/stevor, average var transcription was reduced in 
all transfected lines compared with 3D7. Importantly, 
however, \/ar downregulation was much less pronounced 
in 3D7/phist, 3D7/pfmc-2tm and 3D7/cam compared with 
the lines carrying activated \/ar promoters. Compared with 
3D7/cam, significant downregulation occurred specifically 
only in 3D7/upsA, 3D7/upsB, 3D7/upsC and 3D7/rifA1 
parasites. This is also illustrated by comparing the relative 
expression levels for each individual var gene (Fig. 4B). In 
3D7/upsA, 3D7/upsB and 3D7/rifA1 transcription of all var 
genes was largely reduced. In 3D7/upsC most var tran- 
scripts were also reduced; however, transcripts of the 
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Fig. 4. The capability to interfere with mutually exclusive var transcription is a specific feature of var gene promoters. 

A. Knock-down of endogenous var gene transcription. The heat-map reflects the mean relative expression of all var genes at four time points 
in each cell line (TP1-TP4). The box plot visualizes the distribution of relative var expression levels in all cell lines at the time of peak 
transcription (TP1). Parasites transfected with var promoter constructs are shaded. The box reflects the lower and upper quartiles, 
respectively, and the median is indicated by a horizontal line. Outliers are shown as dots. The solid red and black lines are drawn at the 
median for 3D7/cam and 3D7 respectively. Significant changes in relative var expression compared with the control transfectant 3D7/cam are 
indicated by asterisks (P< 0.01). 

B. The heat-map reflects relative transcript abundance for each var gene at all four time points (TP1-TP4) in all cell lines. Arrowheads identify 
two upsA var genes that were further analysed by qRT-PCR (Fig. S5). 



upsAand upsB/A i/ar subgroups exhibited similar relative 
expression when compared with the 3D7 control. We 
would lil<e to point out though that upsA and upsB/A genes 
are hardly expressed in unselected 3D7 parasites (Frank 
etal., 2007; Peters etal., 2007) and are expected only in 
minute amounts in the 3D7 cDNA reference pool. Conse- 
quently, the hybridization signals in 3D7/upsC most likely 
represent transcription slightly above background rather 
than significant mRNA levels. We addressed this possibil- 
ity by quantifying two upsA transcripts by qRT-PCR and 
confirmed that both genes are hardly transcribed in 3D7/ 
upsC as well as in the 3D7/cam control line (Fig. S5). 
However, we cannot exclude that var transcription is not 
always strictly mutually exclusive in all parasites in the 
3D7/upsC population. Indeed, simultaneous expression 
of two var genes in a single parasite has been docu- 
mented recently (Dzikowski etal., 2007; Brolin etal., 
2009; Joergensen etal., 2010). Importantly, however, and 
in striking contrast to the var and rifA1 promoter transfec- 
tants, transcripts of all var genes were detected in 
the 3D7/stevor, 3D7/phist, 3D7/pfmc-2tm and 3D7/cam 
populations. Compared with 3D7, the relative abundance 
of most transcripts was slightly lower and some were 
transcribed at higher levels in these lines. 

In addition to var genes, other heterochromatic genes 
also showed a trend towards downregulation in trans- 
fected versus untransfected parasites, irrespective of the 
promoter used to drive hdhfr-gfp. However, when compar- 
ing global transcription profiles to the control transfectant 
3D7/cam this effect disappears except for var genes, 
which remain specifically downregulated (Fig. 5). We 
therefore conclude that transcription of heterochromatic 
genes is generally lower in episomally transfected 
parasites. This hypothesis is corroborated by the finding 
that 3D7/stevor, the only cell line where integration of the 
plasmid concatamer into the endogenous locus had 
occurred, was the only cell line displaying wild-type var 
transcript levels. 

Unlike the var promoters, neither of the rifA1, stevor, 
phist and pfmc-2tm promoters caused a specific down- 
regulation in transcription of their cognate endogenous 
gene families (Fig. 6). We observed a trend towards 
reduced transcription of endogenous stevorand pfmc-2tm 
transcription in most transfected lines compared with wild- 
type 3D7, but this was also the case in the 3D7/cam 



control and hence not specific to any of the promoters 
tested. These data clearly demonstrate that transcrip- 
tional control of the rif, stevor, phist and pfmc-2tm family 
promoters does not obey the same strategy employed for 
var genes. These findings are in line with results from 
previous studies (Lavazec etal., 2007; Petter etal., 2007; 
Niang etal., 2009) underscoring the hypothesis that sub- 
telomeric non-var families are not expressed by strict 
mutual exclusion. 

In summary, we show that activation of the three main 
var promoter types, and a head-to-head rifA1 promoter, 
results in silencing of endogenous var genes. Promoters 
of other gene families do not have this competence under- 
scoring the specificity of var promoter-mediated interfer- 
ence with mutually exclusive vartranscription. In contrast, 
variegated expression of other gene families occurs inde- 
pendently of promoter activation. 

Stochastic variation or true regulatory cross-talk? 

Variegated expression of non- var families between indi- 
vidual lines was still obvious, and except for the phist 
families (Fig. 6C) we observed altered average expression 
in some transgenic lines relative to 3D7/cam. For instance, 
in late ring stages (TP2) rif and stevor transcription was 
slightly but significantly altered in several transfected lines 
(Fig. 6A and B). Average pfmc-2tm transcription was sig- 
nificantly reduced in 3D7/pfmc-2tm trophozoites but also in 
3D7/stevor and 3D7/phist, showing that this downregula- 
tion was not specific to the active pfmc-2tm promoter. 
Moreover, in late ring stages average pfmc-2tm transcrip- 
tion was significantly lower in 3D7/upsA and higher in 
3D7/upsB. If some or all of these altered profiles are indeed 
directly linked to the activity of particular promoter types 
and/or shared regulatory factors, or if they are simply due 
to stochastic variation on the population level remains to be 
investigated using more targeted approaches. 

Next, we were interested in asking if any of the recently 
described and uncharacterized hyp families (Sargeant 
etal., 2006) were differentially expressed in our trans- 
fected cell lines. While the temporal transcription profile of 
most hyp genes was similar in all lines, we observed a 
striking upregulation in rtyp4 transcription in all time points 
in 3D7/pfmc-2tm parasites compared with all other lines 
(Fig. S6). This is particularly interesting given that all hyp4 
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Fig. 5. The scatter plots show the correlation of relative expression levels of all genes at TP1. Pearson correlation coefficients are indicated. 
var genes are highlighted in blue, heterochromatic genes in red (Flueck etal., 2009), all other genes are shown in grey. Values on the y- and 
X-axis represent log2 ratios for each gene in the two cell lines. 

A. 3D7/upsB (top), 3D7/pfmc-2tm (centre) and 3D7/cam (bottom) in comparison with 3D7 wild-type parasites. 

B. 3D7/upsB (top) and 3D7/pfmc-2tm (bottom) in comparison with the 3D7/cam control transfectant. 



paralogues are located in head-to-tail orientation directly 
upstream of a pfmc-2tm locus. These observations 
suggest that expression of both families is co-regulated 
and that their exported gene products may cooperate in 
the same functional pathway in the iRBC. As a first 



step towards testing this hypothesis we generated a 
double transgenic cell line simultaneously expressing 
C-terminally tagged PFMC-2TM-GFP and HYP4- 
CherryFP and demonstrate that indeed HYP4 colocalizes 
with PFMC-2TM to MCs (Fig. S6). 
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Lastly, we addressed the interesting question as to 
whether var gene activation involves a single regulatory 
pathway, or if distinct transcriptional regulators mediate 
subtype-specific activation of var genes. To do this, we 
interrogated our datasets to identify genes with differential 
expression in the 3D7/upsA, 3D7/upsB and 3D7/upsC 
transfectants. This comparative approach identified only 
minor transcriptional differences between these cell lines. 
We pinpointed a total of only 20 transcripts with greater 
than threefold changes in relative expression in TP1 and 
TP2 in any one compared with both other cell lines. Eight, 
three and nine genes were differentially expressed 
in 3D7/upsA, 3D7/upsB and 3D7/upsC respectively 
(Fig. S7). Nearly all of these changes affected heterochro- 
matic genes that are expected to undergo variegated 
expression. Interestingly, one of the genes upregulated in 
3D7/upsA encodes the ApiAP2 factor PfSIP2 that has been 
implicated in silencing of subtelomeric upsB-type var 
genes (Flueck etal., 2010). However, apart from this inter- 
esting association this preliminary analysis failed to detect 
a clear link between distinct signatures of potential regula- 
tory factors and activation of different var promoter sub- 
types, suggesting their activation is carried out by a 
common regulatory pathway. 

Discussion 

var upstream regions are the dominant functional ele- 
ments involved in the three main pillars of var gene regu- 
lation - silencing, activation and mutually exclusive locus 
recognition (Calderwood et al., 2003; Voss etal., 2003; 
2006; 2007; Gannoun-Zaki etal., 2005; Dzikowski etal., 
2006; 2007; Frank etal., 2006; Chookajorn etal, 2007). 
Here, we were able to expand our knowledge on this 
complex control strategy onto a genome-wide level, and 
to incorporate the analysis of upsA and several additional 
non-i/ar promoters. Our comparative functional analysis 
separates the i/arand non-var subtelomeric gene families 
into two groups based on a different logic of transcrip- 
tional regulation. 

The three main var promoter types share several func- 
tional attributes including similar relative activities and 
similar temporal activation profiles restricted to ring 
stages. Furthermore, var promoters adopt two different 
states of activity; silenced or active. While this behaviour 
has been described previously for upsB and upsC pro- 
moters, our results suggest for the first time that upsA 
promoters are regulated similarly (Dzikowski etal., 2006; 
Voss etal., 2006; 2007; Muhle etal., 2009). The relatively 
low level of induction observed for the upsA promoter in 
comparison with upsB and upsC has to be interpreted 
carefully. First, for technical reasons upsA-driven tran- 
scripts in unselected (3D7/upsAt) and selected (3D7/ 
upsA) parasites differ in their coding {hdhfr-gfp versus 



bsd) and 3' UTR (hrp2 versus PbDT terminator) 
sequences. These differences may cause differential 
regulation at the level of mRNA stability, which in turn 
would affect the abundance of steady-state transcripts. 
Second, it is impossible to control or assess the propor- 
tion of parasites carrying active promoters prior to selec- 
tion; differences in these proportions are probably 
reflected in the different default activity levels in 3D7/ 
upsAf, 3D7/upsB and 3D7/upsC, which were lowest for 
upsB and highest for upsC. These uncertainties notwith- 
standing, activation of episomal upsA promoters led to a 
specific downregulation of endogenous var transcription 
as previously reported for upsB and upsC. We therefore 
believe that although the active upsA promoter was 
placed in a different plasmid context it still functioned as a 
regulatory entity equivalent to upsB and upsC promoters. 
Hence, despite considerable sequence variation between 
different var promoter types mutually exclusive activation 
of all var genes appears to be governed by a common 
regulatory mechanism. This assumption is also supported 
by our failure to identify distinct signatures of transcripts 
coding for potential transcriptional regulators in 3D7/ 
upsA, 3D7/upsB and 3D7/upsC parasites. 

Our initial attempts to select for an active upsA pro- 
moter by WR selection were unsuccessful. Likewise, a 
previous study failed to obtain activated upsA promoters 
based on drug selection (Muhle et al., 2009). The fact that 
the relative hdhfr-gfp mRNA levels in unselected 3D7/ 
upsAf parasites were higher than those in most other 
WR-selected lines (Fig. S4) demonstrates that the refrac- 
toriness of 3D7/upsAt to WR challenge was not due to 
irreversible silencing. One possibility is that the upsA 5' 
untranslated region (5' UTR) mediates inefficient transla- 
tion, which in turn may result in insufficient hDHFR con- 
centrations to support WR resistance. Interestingly, the 5' 
UTR of var2csa has been implicated in translational 
repression (Mok etal., 2008; Amulic etal., 2009), and we 
recently made similar observations for the upsG 5' UTR 
(unpublished). Alternatively, it is also plausible that the 
phenotypic background in our 3D7 bulk population is 
rather incompatible with expression of upsA genes. In 
unselected cultured lines upsA genes are preferably 
silenced and have faster switch off rates than upsB and 
upsC genes (Dahlback etal., 2007; Frank etal., 2007; 
Peters etal., 2007). It would therefore be interesting to 
test if parasites selected for expression of upsA variants 
provide an environment suitable for activation of the 
pBupsAt construct. Enrichment for upsA-expressing cul- 
tured parasites can be achieved by repeated selection 
with antisera directed against upsA-type PfEIVIPI variants 
(Jensen etal., 2004; Wang etal., 2009). Similarly, upsA 
var subtypes are also frequently expressed in clinical 
isolates and are associated with severe disease (Kirch- 
gatter and Portillo, 2002; Bull etal., 2005; Kyriacou etal.. 
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2006; Rottmann efa/., 2006; Peters etal., 2007). This 
example suggests that environmental cues are important 
in determining the predisposition or capacity to express 
certain classes of var genes. The differential role of silent 
information regulator 2A(PfSIR2A) and PfSIR2B in silenc- 
ing of var gene subgroups and neighbouring upsA var and 



rifA1 genes (Duraisingh etal., 2005; Tonkin etal., 2009; 
IVIerrick etal., 2010) may provide an example of how 
such signals could be translated into changes in 
var expression. Likewise, the observed increase in tran- 
scription of pfsip2 in the 3D7/upsA population may indi- 
cate a role in modulating var transcription. For instance, 
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Fig. 6. Multigene families other than var are not subject to promoter-induced mutual exclusion. Heat-maps reflect the mean relative 
expression of all members of the endogenous gene families at four time points in each cell line. Box plots visualize the distribution of relative 
expression levels for each gene family in all cell lines at the time of peak transcription. The box reflects the lower and upper quartiles, 
respectively, and the median is indicated by a horizontal line. Outliers are shown as dots. The solid red and black lines are drawn at the 
median for 3D7/cam and 3D7 respectively. Parasite lines transfected with the cognate promoter construct are shaded. Significant changes in 
relative expression compared with the control transfectant 3D7/cam are indicated by asterisks (P< 0.01). ns, not significant. 

A. Relative rif expression. 

B. Relative sfevor expression. 

C. Relative phist expression. 

D. Relative pfmc-2tm expression. 



increased nuclear levels of PfSIP2 may be required for 
enhanced silencing of upsB genes in upsA-expressing 
parasites. 

Analysis of the head-to-head rifA1 pronnoter delivered 
several lines of evidence suggesting co-regulation of rifA1 
and upsA var genes. First, the temporal activity profile of 
this rif promoter matches that of endogenous var rather 
than rif genes. Interestingly, a recent study also provided 
evidence for a similar transcription profile (albeit at lower 
temporal resolution) of neighbouring upsA var and rifA1 
genes (Wang etal., 2009). Second, both promoters 
display increased activity levels in selected compared 
with unselected parasites. Third, and most compelling, 
selection for rifA1 promoter activation induced an endog- 
enous var knock-down similar to that achieved by var 
promoters. Overall, these results imply that activation of 
this subset of rif genes may be carried out by the same 
mechanisms and transcription factors that are also 
required for mutually exclusive activation of upsA genes. 
Interestingly, head-to-head upsA and rifA1 genes were 
also co-regulated in parasites selected with IgG antibod- 
ies from highly immune adults (Wang etal., 2009), and 
are similarly upregulated in absence of PfSIR2A (Durais- 
ingh et ai., 2005; Tonkin et ai., 2009). In light of this appar- 
ent co-regulation in vitro and in vivo it will be interesting to 
test if upsA PfEMPI and A1-type RIFINs also cooperate 
on the protein level. If the functional properties of the rifA1 
promoter analysed here is representative for the entire rif 
family remains to be investigated. Based on sequence 
features and subcellular localization patterns, the RIFIN 
family is grouped into several distinct subgroups with 
likely diverse functions, one of which consists of rifA1 
genes positioned head-to-head with upsA var genes 
(Petter etai., 2007; 2008; Joannin etal., 2008; 2011; 
Wang etal., 2009). It is therefore conceivable that func- 
tionally distinct rif gene subsets are controlled by alterna- 
tive regulatory mechanisms. It this context, it has to be 
noted that high copy numbers of plasmids carrying a 
different rrf promoter type also caused downregulation of 
i/ar transcription (Howitt etal., 2009). However, due to 
differences in the experimental set-up of that study com- 
pared with ours it remains unknown whether this effect 
can be attributed to true interference with mutually exclu- 
sive i/ar transcription (see below). 



The stevor, phistb and pfmc-2tm promoters showed no 
evidence for reversible states of activity indicating that 
these promoters are not directly targeted by the silencing 
machinery. Both the stevor and pfmc-2tm promoters were 
properly controlled in terms of temporal activity during the 
IDC. The phistb promoter was active coincident with 
endogenous p/7/sf transcription in ring stages and tropho- 
zoites, but remained active in early schizonts when 
endogenous phist transcription was reduced. If this dis- 
crepancy is due to the absence of critical c/s-acting ele- 
ments on the 2.5 kb episomal promoter, or to differences 
in mRNA stability between the episomal and endogenous 
transcripts is unclear. However, in light of the substantial 
activity displayed by the phistb promoter we envisage it 
was still controlled by the same transcription factors that 
also mediate endogenous phist transcription in ring 
stages and trophozoites. Importantly, for all three promot- 
ers we observed no specific downregulation in transcrip- 
tion of the endogenous rif, stevor, phistb and pfmc-2tm 
families. In these parasites the corresponding plasmid- 
encoded and endogenous promoters were active simul- 
taneously, rather than in a mutually exclusive manner. In 
summary, these results demonstrate that transcriptional 
regulation of non-varfamilies is controlled by mechanisms 
different from those governing mutually exclusive var 
transcription. They furthermore suggest that transcription 
of non-var families is not mutually exclusive, which is 
consistent with previous reports showing that more than 
one RIFIN or STEVOR variant can be expressed concur- 
rently in a single parasite (Lavazec etal., 2007; Petter 
etal., 2007; Niang etal., 2009). 

In a recent study, Howitt and colleagues observed that 
in parasites carrying high copy numbers of unpaired var, 
rif stevor or pfmc-2tm promoters transcription of endog- 
enous var, stevor and pfmc-2tm genes was decreased 
compared with wild-type parasites and to parasites carry- 
ing low copy numbers (Howitt etal., 2009). This down- 
regulation occurred in a family-transcending manner and 
to variable degrees depending on the promoters and gene 
families analysed. To explain these observations, the 
authors proposed the existence of a common titratable 
nuclear factor in activating subtelomeric gene families 
(Howitt etal., 2009). This hypothesis was based on 
several important results obtained by the same group in 



© 2012 Blackwell Publishing Ltd, Molecular Microbiology, 84, 243-259 



254 K. Wifraer et al. ■ 



earlier studies. Episomal var promoters were only 
silenced and counted by the mutual exclusion system if 
paired with a second promoter such as the var intron or a 
heterologous promoter (Deitsch ef a/., 2001; Calderwood 
efa/., 2003; Frank efa/., 2006; Dzikowski efa/., 2007). 
Yet, unpaired varpromoters still caused downregulation of 
endogenous var transcription if present in high plasmid 
copy numbers, but this effect was attributed to absorption 
of a titratable activation factor rather than true infiltration 
of the episome into the mutual exclusion system 
(Dzikowski and Deitsch, 2008). Here, by using an 
improved experimental set-up (including a control trans- 
fectant, paired promoters, sampling of consecutive IDC 
time points and genome-wide analyses) we are now able 
to revisit the above hypothesis and to distinguish between 
specific and non-specific effects. Similarly to Howitt and 
colleagues, we observed reduced transcription of var 
genes and other subtelomeric gene families in all trans- 
fected cell lines. However, the only specific effect 
occurred in 3D7/upsA, 3D7/upsB, 3D7/upsC and 3D7/ 
rifA1 where paired var and the r/M7 promoters caused 
specific silencing of endogenous vartranscription by inter- 
ference with mutually exclusive transcription. In contrast, 
the moderate degree of var downregulation in the remain- 
ing transfectants, and the reduced transcription of stevor 
and pfmc-2tm genes, is a likely result of non-specific 
mechanisms related or identical to those observed using 
unpaired promoter constructs (Dzikowski and Deitsch, 
2008; Howitt efa/., 2009). Notably, this effect also 
occurred to the same extent in the control transfectant 
3D7/cam. Hence, our results are inconsistent with absor- 
bance of a limiting factor specifically by subtelomeric 
virulence gene promoters; they rather indicate that subte- 
lomeric gene transcription is generally reduced in trans- 
fected parasites. At this stage we can only speculate 
about the nature of this non-specific effect. It is conceiv- 
able that high plasmid copy numbers would absorb limit- 
ing factors critical for the expression of heterochromatic 
genes such as histone variants, chromatin remodelers 
and/or general transcription factors. Alternatively, reduced 
transcription of non-essential genes in transfected popu- 
lations may be a direct or indirect consequence of mecha- 
nisms that compensate for a possible reduction in 
parasite fitness due to drug selection and forced mainte- 
nance of large concatenated plasmids. The transcriptional 
changes commonly observed in all transfected lines may 
be a reflection of such compensatory mechanisms. 

A currently accepted model proposes the existence of a 
unique and physically restricted perinuclear zone dedi- 
cated to the expression of a single var locus only (Durais- 
ingh efa/., 2005; Ralph efa/., 2005; Voss efa/., 2006; 
Dzikowski efa/., 2007; Lopez-Rubio efa/., 2009). Accord- 
ing to this concept, transcriptional var switching would 
occur through competition between a silenced and the 



active var gene. Our results are concordant with this 
hypothesis. Among the gene families investigated here, 
the mechanism by which singular var gene choice occurs 
appears to be unique and specific to the var gene family; 
promoters of other gene families do not interfere with this 
system. We envisage that specific c/s-acting elements in 
var upstream regions could target var loci to such an 
expression site and that default silencing may be a pre- 
requisite to protect var loci from illegitimate expression 
elsewhere in the nucleus. In fact, the reversible respon- 
siveness of var promoters to silencing and activation may 
represent an inherent contributor to the overall mecha- 
nism controlling mutually exclusive transcription. On the 
contrary, r/f, stevor, phistb and pfmc-2tm promoter activity 
is not mutually exclusive and appears to be mediated by 
an existing pool of readily available specific transcription 
factors. Recent genome-wide ChIP experiments demon- 
strated the ubiquitous enrichment of H3K9me3 and 
PfHP1 over P. fa/c/parum-specific gene family islands 
throughout the genome (Flueck efa/., 2009; Lopez-Rubio 
efa/., 2009; Salcedo-Amaya efa/., 2009). var gene acti- 
vation is linked to the local exchange of H3K9me3/PfHP1 
with H3K9ac and H3K4me2/3 predominantly in the 
upstream region (Lopez-Rubio efa/., 2007). Such a func- 
tional epigenetic footprint is probably also involved in tran- 
scriptional activation of silenced non- var family members. 
However, the apparent lack of autonomous c/s-acting ele- 
ments capable of recruiting the silencing machinery to the 
stevor, pfmc-2tm and phistb promoters suggests that their 
expression may be passively controlled by cis and/or 
trans effects of neighbouring sequences that actively par- 
ticipate in structuring the local chromatin environment. 

Overall, our study demonstrates a lack of systematic 
regulatory cross-talk between different subtelomeric gene 
families in P. falciparum. The most interesting circumstan- 
tial evidence for potential co-regulation is related to the 
pfmc-2tm family. First, we observed a significant down- 
and upregulation of pfmc-2tm genes in 3D7/upsA and 
3D7/upsB respectively. It is worth mentioning that similar 
to our in vitro data a previous study reported inordinate 
pfmc-2tm expression patterns in vivo (Mackinnon efa/., 
2009). However, if pfmc-2tm expression is truly 
co-regulated with subtype-specific var gene transcription, 
or if these variations result from stochastic fluctuations in 
temporal pfmc-2tm activity, needs to be tested in future 
studies. Second, iiyp4 genes displayed a striking upregu- 
lation in 3D7/pfmc-2tm compared with all other cell lines. 
This result, together with the strict positional link between 
both gene types in the genome, and the colocalization of 
PFMC-2TM and HYP4 to MCs suggests that these two 
gene families may be functionally inter-connected. 

To our knowledge, this is the first study comparing 
multiple transfected P. falciparum lines by global tran- 
scriptional profiling. As an important notion for future 
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Studies, we advise to use mocl< transfectants ratlier than 
untransfected parasites as controls for comparative tran- 
scriptome analyses involving transgenic lines. Another 
important technical aspect relates to the use of transgenic 
P. falciparum to study the function of variant surface 
antigens. While the ability to create PfEMPI knock-down 
parasites using such a transfection-based approach con- 
tinues to be informative for the analysis of immune 
responses to and binding properties of PfEMPI (Voss 
etal., 2006; D'Ombrain etal., 2007; Elliott etal., 2007; 
Maier etal., 2008), this approach may not be feasible to 
achieve a complete knock-down of other suspected 
surface antigens and virulence factors. 

In this work we shed important light on the transcrip- 
tional control of virulence gene families in P. falciparum 
and possibly that of other pathogens. We show that acti- 
vation of the three major types of var gene promoters is 
equally important in the process of singular var gene 
choice in P. falciparum. The transition of a var gene pro- 
moter from the silenced to the activated state most likely 
represents an integral part of the overall mutual exclusion 
mechanism. We also demonstrate that this mode of regu- 
lation is unique to var genes; non-\/ar families such as rif 
stevor, phist and pfmc-2tm do not obey the same principle 
of mutually exclusive transcription. This difference in the 
control of i/ar versus non-vargene families underscores 
the high significance of antigenic variation of PfEMPI in 
parasite survival and transmission. 

Experimental procedures 

Parasite culture and transfection 

Plasmodium falciparum 3D7 parasites were cultured as 
described previously (Trager and Jenson, 1978). Transfec- 
tions were performed as described (Voss etal., 2006) and 
parasites were selected on 2.5 |xg ml"' BSD and 4 nM 
WR99210. Generation of transfection constructs is explained 
in detail in Supporting Information. Growth synchronization 
was achieved by repeated sorbitol lysis (Lambros and 
Vanderberg, 1979). To obtain an 8 h growth window pre- 
synchronized parasites were synchronized twice 16 h apart. 

Transcriptional profiling by microarray analysis 

For microarray analysis, total RNA was harvested at four 
consecutive time points: 6-14 hpi, 14-22 hpi, 22-30 hpi and 
30-38 hpi. A pool of total RNA isolated from synchronized 
3D7 parasites at five consecutive time points (10 h growth 
window each) covering the complete intra-erythrocytic devel- 
opmental cycle was used as reference sample. Total RNA 
was prepared by disrupting 1 .5 ml of RBC pellet (> 3% para- 
sitaemia) with 10 ml of TriReagent (Sigma-Aldrich) following 
manufacturer's instructions. cDNA was generated from total 
RNA, and the test and reference samples were labelled as 
previously described (Bozdech etal., 2003) and hybridized to 
a long oligonucleotide genome-wide microarray containing 



10 166 probes (Hu etal., 2007). The raw microarray data 
representing mRNA abundance ratios between each sample 
and the reference pool were subjected to lowess slide nor- 
malization and background filtering as implemented by the 
Acuity 4.0 program (Molecular Devices). The Cy5/Cy3 
log ratios of multiple probes per gene were averaged. Genes 
with detectable signals in less than 75% of all samples as well 
as tRNA and rRNA genes were excluded from the analysis. 
All microarray data are available at the GEO repository 
(Accession No. GSE31829). 

Data analysis 

Genes and arrays were clustered by average linkage clus- 
tering using Gene Cluster 3.0 (Eisen etal., 1998). The array 
tree and heatmaps were generated using the Java Treeview 
program (Saldanha, 2004). The cut-off for changes in relative 
expression was fourfold for pairwise comparisons (Fig. 2), 
and threefold (Fig. S7) or twofold (Fig. S2) to identify genes 
differentially expressed in one versus several other cell lines. 
Boxplots and scatter plots were generated using the statisti- 
cal package R (version 2.10.0). Paired and unpaired Mests 
were used to assess statistical significance between means 
of log2 ratios obtained from microarray hybridizations and 
qRT-PCR triplicate values respectively. 

Genomic DNA extraction, RNA extraction and cDNA 
synthesis 

Parasites were released from iRBCs in 0.15% saponin. Total 
RNA was isolated using TriReagent (Sigma-Aldrich) and 
further purified using the RNeasy Plus Mini Kit (Qiagen) for 
removal of gDNA. Residual gDNA was digested with TURBO 
DNA-free™ (Ambion). All samples were tested negative for 
contaminating gDNA by quantitative PGR (qPGR). RNA was 
reverse transcribed using the RETROscript Kit® (Ambion). 
gDNA was extracted from the same parasite samples by 
addition of 3% SDS/0.27 M Na-acetate to the parasite pellet 
followed by phenol/chloroform extraction and precipitation. 
gDNA was resuspended in ddH20 and diluted 1:100 for 
qPCR analysis. 

Quantitative real-time PCR 

For determination of the temporal promoter activity profiles 
RNA and gDNA were harvested from the same cultures that 
were also used for microarray analysis. For analysis of pro- 
moter activities in unselected and drug-selected populations, 
and to assess upsA transcripts, parasites were synchronized 
as described above and gDNA and total RNA was harvested 
at 6-14 hpi (TP1) for 3D7/upsA, 3D7/upsB, 3D7/upsC, 3D7/ 
rifAl, 3D7/phist and 3D7/cam, and at 22-30 hpi (TP3) for 
3D7/stevor and 3D7/pfmc-2tm respectively. qRT-PCR experi- 
ments were performed on three independent biological 
samples. Plasmid copy numbers were determined by qPCR 
on gDNA and calculated by dividing the absolute copy 
numbers obtained for hdhfr-gfp by the value obtained for the 
PFL0900C (arginyl-tRNA synthetase) locus (Frank etal., 
2006). All reactions were run in triplicate yielding virtually 
identical Ct values. Five serial 1:10 dilutions of 3D7 genomic 
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or plasmid DNA were used as standard for absolute 
quantification. Relative transcript profiles were calculated by 
normalization against the housekeeping gene PFL0900c. 
Cycling conditions were: 50°C, 2 min; 95°C, 10 min, followed 
by 40 cycles of 95°C, 15 s; 58°C, 1 min. Product-specific 
amplification was ensured by performing melting curves for 
each reaction. Reactions were performed at final primer con- 
centrations of 0.4 |iM using Power SYBR® Green Master Mix 
(Applied Biosystems) on a StepOnePlus™ Real-Time PGR 
System (Applied Biosystems) in a reaction volume of 12 nl. 
Primers are listed in Table SI. 

Live cell fluorescence microscopy 

Five hundred microlitres of culture (5% haematocrit, 3-8% 
trophozoites) was pelleted for 1 min at 200 g and washed 
twice in 1 ml of RPMI-HEPES. The pellet was resuspended in 
400 \i\ of RPMI-HEPES containing 2.5 |xg mM DAPI and 
incubated for 10 min at 37°G. Three microlitres of packed 
RBGs were mixed with 2 \i\ of Vectashield (Vector Laborato- 
ries) and immediately analysed on a Leica DM 5000B 
microscope. Images were taken at lOOx magnification with a 
Leica DFG 300 FX camera and acquired via the Leica IM 
1000 software, and processed using Adobe Photoshop GS2. 

Souttiern blotting 

gDNA was subjected to Southern analysis using standard 
procedures. gDNA was digested with EcoRI/Pvull and sepa- 
rated on 0.7% agarose gels. Southern blots were probed with 
random-primed ^^P-dATP-labelled hdhfr. 
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